WELLPOSEDNESS OF THE NAVIER-STOKES-MAXWELL 

EQUATIONS 



PIERRE GERMAIN, SLIM IBRAHIM, AND NADER MASMOUDI 



Abstract. We study the local and global wellposedness of a full system of 
Magneto-Hydro-Dynamic equations. The system is a coupling of the forced 
(Lorentz force) incompressible Navier-Stokes equations with the Maxwell equa- 
tions through Ohm's law for the current. We show the local existence of mild 
solutions for arbitrarily large data in a space similar to the scale invariant spaces 
classically used for Navier-Stokes. These solutions are global if the initial data 
are small enough. Our results not only simplify and unify the proofs for the space 
dimensions two and three but also refine those in [TO]. The main simplification 
comes from an a priori Lf estimate for solutions of the forced Navier-Stokes 
equations. 



1. Introduction 

The purpose of this paper is the study of the following full Magneto-Hydro- 
Dynamics system (MHD): 

f^r + v ■ Vv - uAv + Vp = jxB 

d t E — curlB = —j 

(1.1) < d t B + cur\E = 

divv = divB = 

a(E + vxB) = j 

subject to the initial data 

v lt=0 = v°, B {t=0 = B°, E {t=0 = E°. 

Here, v,E,B : Rf x R% — > R 3 are vector fields defined on R d (d = 2 or 3). The 
vector field v = (v±, ...,t>d) is the velocity of the fluid, v its viscosity and the scalar 
function p stands for the pressure. The vector fields E and B are the electric and 
magnetic fields, respectively, and j is the electric current given by Ohm's law (the 
fifth equation of the system, where a is the electric resistivity). The force term 
j x B in the Navier-Stokes equations comes from the Lorentz force under a quasi- 
neutrality assumption of the net charge carried by the fluid. Note that the pressure 
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p can be recovered from v and j xB via an explicit Calderon-Zygmund operator (see 
[5] for instance). The second equation in (jl.lj) is the Ampere-Maxwell equation for 
an electric field E. The third equation is nothing but Faraday's law. For a detailed 
introduction to MHD, we refer to Davidson [8] and Biskamp [2]. 

Note that in the 2D case, the functions v, E, B, and j are defined on the whole 
space M. 2 with values in M 3 . In this case, the operator V is given by 

V = (d Xl ,d X2 ,0) T . 

Thus 

diw := d Xl Vi + d X2 v 2 , Vp := (d Xl p, d X2 p, 0) T , 

and 

curlF := {d X2 F 3 , -8 Xl F 3 ,d Xl F 2 - d^F x ) T . 



In the following, we take a = v = 1 to alleviate the notations. 

Multiply the Navier-Stokes equations in (jl.ip by v, the Ampere-Maxwell equations 
by (B,E) T and integrate (using the divergence free condition on the velocity); this 
gives the formal energy identity 

^[11* + W B Wl> + W E \\h] + WiWh + HV«||i a = o. 

This identity shows that the energy is dissipated by the viscosity and the electric 
resistivity. It also suggests that one should be able to construct a global finite energy 
weak solution (a la Leray) for data lying in L 2 (IR rf ). However, this intuitive expec- 
tation remains an interesting open problem for in both dimensions d = 2,3. 
Indeed, given a standard approximating scheme, it is hard to obtain compactness of 
the solutions, especially for the magnetic field due to the hyperbolicity of Maxwell's 
equations. In dimension 2, the equation is energy critical, but running a fixed point 
argument for data (v°, E°, B°) only in L 2 (W d ) 3 seems very difficult due to the term 
Ex B. 



Imposing more regularity on the initial electro-magnetic field, existence results 
are known. Recently, for initial data (i; ,^ ,^ ) G L 2 (R 2 ) x (H S (R 2 )) 2 with s > 0, 
Masmoudi in [JJ] proved the existence and uniqueness of global strong solutions to 
(jl.ip . His proof relies on the use of the energy inequality combined with a logarithmic 
inequality that enabled him to upper estimate the L°° norm of the velocity field by 
the energy norm and higher Sobolev norms. It is also interesting to note that the 
proof in [14] does not use the divergence free condition of the magnetic field, nor 
the decay property of the linear part coming from Maxwell's equations, namely 



(1.2) 



dE 

at 



curl B + E 



H + curl E 



f. 



0. 



V 



V-B 



= 0. 
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Another line of research was pursued by Ibrahim and Keraani [TO] who considered 
data (v ,E ,B ) € B^(R 3 ) x (ij 1 / 2 ^ 3 )) 2 in dimension d = 3, and (v ,E ,B ) e 

#2,1 ( r2 ) x (- l2 o 9 ( r2 )) in dimension (f = 2 (see below for the definiton of these 

functional spaces). These authors built up strong solutions by using parabolic reg- 

ularization arguments giving control of the L°° norm of the velocity field of the 

solution. More recently, Ibrahim and Yoneda constructed local in time solution for 

non-decaying initial data on the torus. See [11] for more details. 

In this paper, we follow up on the work of Ibrahim and Keraani by running a 

fixed point argument to obtain mild solutions, but taking the initial velocity field 

-—1 

in the natural Navier-Stokes space H2- 1 . Our main theorem extends the earlier 
results that were mentioned in many respects: the regularity of the initial velocity 
and electromagnetic fields is lowered, and we unify the proofs in the cases of space 
dimension 2 and 3. One of the key ingredients will be to use an L 2 L°° estimate on 
the velocity field, which simplifies greatly the fixed point argument; in particular, 
the weak decay for the electromagnetic field is not needed any more in dimension 3. 
Before stating our main result, we need a few definitions. 

Definition 1.1. First, let V denote the Leray projection on divergence-free vector 
fields. 

A function T := (v, E, B) with div(v) = div(i?) = is said to be a mild solution 

on a time interval [0, T] of the full MED problem flU} if V € C([0, T], iH" 1 ) and 
satisfies the integral equation 

r(t) = e tA r(0) + f e^-^ A M{T{t')) dt', 
Jo 

with 

I A \ 

^= —I curl 
\ -curl / 

and M{T) = [V(-V{v ® v) + E x B + (v x B) x B), -v x B, 0) T . 

The functional analytic framework we will use is the following. 

Definition 1.2. Let A q denote the dyadic frequency localization operator defined in 
section 2. For s,f£K and a > define the space H a ' by its norm 

q<0 q>0 

We will also use the short-hands 

H S = H S Q ' S , Hf og :=H s { s and H s ' l :=H^. 

Finally, define LZpHa 1 by its norm 

with obvious generalizations to U T H S etc... 
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The space Hf is articulated on the standard homogeneous Sobolev space H s 

with an extra logarithmic weight for the high frequency part. The space H S:t is 
nothing but the usual Sobolev space H* for high frequencies while it behaves like 
H s for low frequencies. If s > t, it is not difficult to see that H s ' = H s + H t . The 
L spaces were first used by Chemin and Lerner [7]. 
Our main result can be stated as follows. 

Theorem 1.3. • In dimension two and for any 

T° := (v°,E°,B°) G L 2 (R 2 ) x L 2 og (R 2 ) x L 2 og (R 2 ), 

there exists T > and a unique mild solution T = (v, E, B) of with 
initial data T° and 

v G L°°(0, T; L 2 ) n L 2 (0, T; H 1 n L°°) 

E G L°°(0, T; Ly og ) n L 2 (0, T; L 2 og ) 

5 G Z oo (0 ) T;Lg )g )nL 2 (0,r;fl' 1 ' ). 

Moreover, the solution is global (i.e. T = oo) if the initial data is sufficiently 
small in L 2 x L 2 x L 2 . 
• in dimension three and for any 

T° := (v°,E°,B°) G i?5(R 3 ) x if^(M 3 ) x if^(]R 3 ) 

t/iere exists T > and a unique mild solution T = (y, E, B) of (jl.ip mi/i 
initial data T° and 

v G L°°(0,r;^l) nL 2 (0,T;ijt nL°°) 

£ G L°°(0,T;H^) nL 2 (0,T;ij^) 

5 G L°°(0,T;iH) nL 2 (0,T;i^l'5). 

Moreover, the solution is global (i.e. T = oo) if the initial data is sufficiently 
small in x x H% . 

In dimension two, the extra logarithmic regularity is needed to estimate the term 
E x B appearing in the Navier-Stokes equations. 

In dimension three, the control of B in L 2 (0, T; H 2>2 ) is not needed to close the 
fixed point estimate, but we added it for completeness. 

The rest of this paper is organized as follows. In the next section we define 
some further tools needed in the proof. In Section 3, we detail the linear (parabolic 
regularity) and nonlinear (product law) estimates needed in the proof of the main 
theorem. The main theorem is then proved in Section 4. Finally, the proofs of some 
technical estimates are given in the appendix. 
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2. Notations and functional spaces 

Throughout this work we use the following notation. 

(1) For any positive A and B the notation A < B means that there exists a 
positive constant C such that A < CB. 

(2) c will always denote an absolute constant < c < 1. 

(3) For any tempered distribution u, both u and Tu denote the Fourier transform 
of u. 

(4) For every pG [1, oo], || • denotes the norm in the Lebesgue space LP . 

(5) For any normed space X , the mixed space-time Lebesgue space L p ([0, T], X) 
denotes the space of functions / such that for almost all t G (0, T), f(t) G X 
and \\f(t)\\ x G LP{Q,T). The notation D>([0,T\,X) is often shortened to 

LlrpX . 

Let us recall the well-known Littlewood-Paley decomposition and the correspond- 
ing cut-off operators. There exists a radial positive function ip G D(K d \{0}) such 
that 



q& 



Supp c^(2~ q -) n Supp tp(2- 



V(eR rf \{0}, 
Vlff-il > 2. 



For every g G Z and i; G 5 



//W\ 



we set 



j=-oo 

Bony's decomposition [3] consists in splitting the product uv into three parts^: 

uv = T u v + T v u + R(u, v), 



with 



T u v = ^Sg^uAgV, R(u,v) = ^AguAyV and A q = ^ A q+i . 



For (p, r) G [1, +oo] 2 and s G M we define the homogeneous Besov space Bp r as the 
set of u G <S'(R d ) such that u = ^ q A q u and 



(2^||A Q n|| iP ; 



e r (z) 



< oo. 



In the case p = r = 2, the space -B| 2 turns out to be the classical homogeneous 
Sobolev space H s . Finally, define L3pBp r is given by distributions u such that 



\U\ rq 



T p,r 



^nU\\ T q TP 



e r (z) 



< 00. 



1 It should be said that this decomposition is true in the class of distributions for which 
^ ?gz A 9 = /. For example, polynomial functions do not belong to this class. 
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3. Linear and Nonlinear estimates 

We will make an extensive use of Bernstein's inequalities (see [5] for instance). 

Lemma 3.1 (Bernstein's lemma). There exists a constant C such that for any 
q,keN,l<a<bandforfe L a {R d ), 

sup \\ff*S q f\\ L > < C k 2^ k+d ^^\\S q f\\ L «, 

\a\=k 

C- k 2i k \\& q f\\ L a < sup\\d a A q f\\ La <C k 2" k \\A g f\\ L a. 

\a\=k 

The parabolic regularity result we will need reads 

Lemma 3.2 (Parabolic regularization, see for example [TJ). Let u be a smooth 
divergence free vector field solving 

dtu — An + Vp = f 

U\ t=0 = u°, 

on some time interval [0, T]. Then, for every p > r > 1 and s € K and j > 1, 

ll U ll s+ 2 < \\u°\\r,s +II/L . s _9+2- 

The following result is an L^L°° estimate which was originally proved in [T2], [B] 
in dimension two. 

Lemma 3.3 (L 2 L°° estimate). Let d = 2,3 andu be a smooth divergence free vector 
field solving 

' d t u -Au + Vp = f 1 + f 2 
«|t=o = u °, 

on some time interval [0,T]. Assume that f\ £ L\,Hz~ l and fi G L^B^i ■ Then, 



(3-1) \H\l%l°° % IK||^d_i + \\h\\ Llk d-i + 11/2 

Proof. Thanks to Lemma 13.21 and using the embeddings 

L 2 Bj ^ L 2 B$? ^ L 2 L°°, 
we can assume that fi = 0. Duhamel's formula gives 

u (t) = e tA u + [ t e^ A Vf 1 (t')dt r , 
Jo 

and thus 

fT 

(3.2) \\u(t)\y T L°° < \\e tA u \\ L ^ + J \\e tA Vh(t')\\ L 2 TLao dt' . 

Using the embedding ija" 1 <^4 B^ 2 an d the characterization of Besov spaces of 

negative regularity (see for example [I]), 

II II 11 11 tA n n 

" "B«2 ~ U H L °°IIl 2 (0,oo)» 



thus we obtain ()3.1|) as desired. □ 
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Now we focus on Maxwell's equations. The first result is an energy type estimate. 



d 

2 ' 



Lemma 3.4. Let a > 0, G\ € L^H^ , and (E,B) be a smooth solution of 

d t E - curl B + E = G 
d t B + curl E = 

(E, B)\ t=Q = (E ,B ), 

on some time interval [0,T]. Then, the following estimate holds (with constants 
independent ofT) 

(3.3) \\E\\_ dj 9> 4_i + ||B||. . d_! < ||(£ ,-Bo)|| . d_ 1 + ||G|| d_ 1 . 

T oo tt 2 ,—, r 2 tt 2 J oo tt 2 tt 2 t 2 tt 2 

l-Jrp Oq I I Lj rp 1L q, I_Jr-p LL Q, ii Q, lj rjp JTL (y 

Moreover, B satisfies the following decay estimate 

(3-4) ll B ll llft f4-.SII(^,Bo)ll ft| _, + l|G|| i2r . | _ I . 

We emphasize that for the existence and uniqueness part of Theorem 11.31 i n 
dimension 3, estimate (I3.4l)is irrelevant. 



Proof. Only the estimate of ||£|| , d d_ 1 requires a proof, which is given in the Ap- 

L H a 

pendix. All other estimates can be derived by a standard energy estimate; apply A q 
to the system, derive an energy inequality, multiply both members of that inequality 
by 2 9 (^~ 1 ) v / max(l,g a ) and take the ^ 2 (Z)-norm. □ 

The following is a series of nonlinear estimates needed for the contraction argu- 
ment. 

Proposition 3.5. For all smooth functions u, E and B defined on some interval 
[0,T], we have the following estimates, with constants independent of T : in space 
dimension 2, 

(3.5) || V(u®u) ~ II ^ II 2,2 (jt oo (K2)niyl (R2)) 1 1^ 1 1 i|, (K 2 )!"!^ 1 (K 2 )) 

(3.6) \\E X -B||^e 2 ;^R 2 )+Li,L2(R2) < \\ E \\LlLf og (K?)\\ B \\l^Lf og (M 2 )nL^H 1 '0 

(3.7) \\U X B|| L ^ L 2 (R2) < |h||i2( it x.( M 2) ni jl( R 2)) ll-BH^oo^ 



and in space dimension 3, 



(3.8) ||V(u<g> v )ll L i^i (R3) ~ H'"llL2,(L-'(R3)niji(R3))ll t; llL2 i(L o ( R 3)ni/l(R3)) 

(3.9) \\E x £11 i < ||£|| , . i , ||£|L . i , 

(3.10) ||u X ■S||_ L a_Hr^( K 3) ~ " U "L2,(L°°(R3)ni/i(R3))l' jB "L5?i : f^(R3)' 
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Estimates (|3.5f) and (|3.8p enable us to control the advection term in the Navier- 
Stokes component of the system in dimension two and three, respectively. Estimates 
(|3.6j) and ()3. 9f) are needed to control the Maxwell part in the Navier-Stokes compo- 
nent. To estimate the term (ux B) x B, we use fl3.6|) . (13. 7j) in two space dimensions 
and (|3.9p . (13.10p in three space dimensions. 

Remark 3.6. Ignoring the time variable, estimate (13. 9h is a particular case of the 
product law 

H Sl (R d )-H S2 (R d ) ^ B%f"~*(R d ) 

withsi,S2 G] — |, |[ and si + S2 > 0. Indeed, it corresponds to the admissible choice 
si = S2 = 2- However, this product law becomes critical in two space dimensions. 
Estimate (13.6|) shows that an extra logarithmic loss is needed in this case. 

We give the proof of the above proposition in the Appendix. 

4. Proof of Theorem 11.31 

4.1. Small data, global existence. Let a = 1 if d = 2 and a = if d = 3. Let 2 

be the set of F := (it, .E, i?) 7 " such that 

u£Z u := L 2 (0, oo, Hi n n L°°(0, oo, 

£€2 £ := n L 2 )(0, oo, J" 1 ) 

B eZ B := L°°(0, oo, .ffj" 1 ) n L 2 (0, oo, Hi' 1 ' 1 ). 

Endow Z, Z u , Z E and Z B with the natural norms. Recall that we seek a solution 
to (II. ip in the integral form 

r(t) = e M r(o) + / ^'-''Ka/W)) < 

with 

/AO \ 
,4= —I curl . 
\ -curl / 

and J\T(T) = (7>(-V(u 8ti) + £xB|(uxB)xB) r MX 5, 0) T . Let B s be the 
ball of the space Z^ centered at zero and with radius 5 > to be chosen. Define 
the map $ on that ball as follows 

$ : B s C Z — > Z 

(4.i) r h. $(r) := / V-^a/V'" 4 ^ + r(0) dt'. 

Jo 

Claim 4.1. 7/||r°|| d d_ 1 d_ 1 < n5, with 5 > and k > sufficiently small, 
then the map & is a contraction on Bg. 
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The theorem follows immediately from the claim: Picard's theorem gives the 
existence of a fixed point of call it T. Then e tA T° + is the desired solution. 

PROOF of the claim: First, notice that <3? (— e tA T ) = 0, while by Lemmas 1331 13.31 
and 13.4 



(4.2) \\e tA T \\ z <C\\T \\. d . d < Ck5 < S - 



for k small enough. On the other hand, we will prove below that, if T\ and r 2 belong 
to B s , 

(4.3) ||$(r 1 )-$(r 2 )|| z <i||r 1 -r 2 || 2 

under the assumptions of the claim. 

The estimates (|4.2p and (|4.3I) easily yield the claim. 

To prove (|Q|h let := {u^E^Bj) 7, € £5 for j = 1,2. Write further 

e'^ + T^t) = (uj,Ej,Bj) 

and set r := T x - T 2 := (u,E,B) T , and $(Fj) := P,- = (u j ,E j ,B j ) T be given by 
P~T]) . Let f := fi - f 2 := (u,E,B) T . We decompose u into tt = m ws + tt M , with 
u NS accounting for the convection term 

ct 



u NS 



I e ( i - i ')^pv(n ® u x + u 2 ® u) dt', 
Jo 



and u M for the Lorentz force 

ct 



u M 



[ e^ A V{ExB 1 +E 2 xB) dt' 
Jo 

+ J e {t - tl)A v({u x Si) xBi + [i*2 xB]xB! + [u 2 x B 2 ] x B^j dt' . 

Moreover, the electro-magnetic field (E, B) satisfies 

(4.4) d t E - curl B + E = u x B x + u 2 x B 

(4.5) d t B + cuv\E = 

with data. First, by Lemma 13.21 Lemma 13.31 and the embedding L^H 2 -1 
L H2 , we have 

',NS\\ < 



II" \\z* < ||7'V(«®ui + iZ2®tt)|l Ll ^_ 1 
and 

||n M |b«< 11^(^x5! +^2X5)11 # _ 2 „ 

+ ||P ((it x Si) xBi + [u 2 xB]xB l + [u 2 x B 2 ]x B 
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Second, applying estimates (|3.5I) and (|3.8j) . we obtain for the convection term 

\\u N \\zu < \\u\\ . d \\Uj\\ ■ d 

5=1,2 

(4-6) < \\T\\ Z Y, (llr.b + He^rOb), 

5=1,2 

whereas the Lorentz force term can be estimated by (|3,6p . (|3.7p . (|3.9p . and (|3.10p 

(4.7) 

+ ||«xB 1 || d.JI-BiH. .d^ d d ! + ||«2 xB|| dJlJSiH. .d^ „ d 1 

~ ll-^llz^ ll-^i \\z B + ll-^Hzs ll-BH^s + ||«||^||5i|| z b||5i|| z b + ||u2|U«||-B||^s||i?i|| Z B 
+ ll^ll^ll^lbsll^lb-B 

■ ll e 1 Hz "I" \\*-j\\z UMl'^ 

5=1,2 



<l|r|U£ [l|e^r |U + || e ^r°||| + 11^11^ + 11^1 



It remains to estimate the electro-magnetic field components of V. Applying the 
energy and the decay estimates (|3.3p to the system (|4.4p . we get 

2 

(4.8) + ||B|| 2fl < unu £(l|e u r°|| 2 + lir.lb). 

5=1 

Gathering the estimates (|4.6p . (|4.7p and (|4.8p gives 

(4-9) ||f \\ z < \\T\\z (\\e tA T°\\ z + ||e M r°||| + \\T 3 \\ Z + \\T 3 \\l) . 

Choosing 5 small enough gives (|4.3p . 



4.2. The local existence. Decompose the initial data (u°, E°, B°) = (u®,E®, B° r ) + 
(u®,E®,Bg) where (u°, E®,, B®.) is regular (say in H 2 ) and (u° s ,E®,B®) is small in 

■ d ■ — — 1 ■ — — 1 

H2~ l xH% xH£ (this can be done by a Fourier cut-off). We look for a solution 
T of {EEI of the form (it, B) = (it s , £ s , v3 s ) + (u r ,E r , B r ) with 



(4.10) 



subject to the initial data 

U r \t=0 = u r> Br\t=o = B r , E r \ t=0 = E®.. 



^ + U r ■ Vll r - Au r + Vp r 


— jr x Br 


d%E r — curl B r 


= ~jr 


dtB r + curl E r 


= 


divit r = divv3 r 


= 


(E r + u r xB r ) 


= 3r 
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Arguing as in [IT], we know that (|4.10p has a unique regular solution. Now we solve 

for (u s , E s , B s ). We have 

(4.11) 

+ U s • Vli s — Au s + U s ■ Vli r + U r ■ Vli s + S/p s = j X B — j r X B r 

d t E s - curl B s = j - j r 

d t B s + cmlE s = 

divtis = divB s = 

subject to the initial data 

u s \ t=0 = u s , B s \ t=0 = E\ t=Q = Eg. 

Proceeding similarly as for the small data result, set 

. = ^ e (t-t')A ( V [-u s ■ Vu s - u s ■ Vu r - u r ■ Vu s + j x B — j r x B r ] \ ^ 
Jo \ 3 ' ~ Jr J 

where T is defined by 

(u s ,Bs,Es) = e tA (u° s ,B° s ,E° s ) + T. 

Applying the same proof as for the small data existence, we can show that the 
map $ is a contraction if we choose a time of existence T sufficiently small. The 
main difference is that new linear terms (in T) appear in <£. These linear terms need 
to be small (as linear maps) for $ to be a contraction; this can be achieved by using 
the smoothness of B and by choosing T small enough. 

For instance: 
(4.12) 

°* e (t-t>)A VEs xBrdt , 



S \\E. q xB r \\_, .di S, \\E S I _ - d 1 \\B r 



The key point is of course that ||-B r || £ i H 2 can be made arbitrarily small by choosing 
T small enough. 

5. Appendix: Proof of the decay for B and Proposition 13.51 

Recall the main part of Lemma 13.31 We emphasize on the fact that this property 
is an extra information about a weak decay that the magnteic field satisfies and it 
has no impact on the well posedness result. 

■ --1 

Lemma 5.1. Let a > and G 6 L^Ha , and (E, B) be a solution of 

d t E - curl B + E = G, 
d t B + curl E = 

on some time interval [0, T\. Then, the following estimate hold with constants not 
depending upon time 



(5-1) aa, < \\(E ,B Q )\\ . ^ + ||G?| 
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Proof. Because of the divergence free property of B, we have 



(5.2) dttB — AB + d t B = curl G, (B, d t B) ]t=0 = (B°, B 1 ) 

Thus, the magnetic field B satisfies an inhomogeneous damped wave equation 
(|5.2j) . In the sequel we denote by C\{t) and £2^) the propagators associated to the 
Fourier multiplier functions 



*x(t,f) =e-/ 2 cosh ( yT74^Ft), < Ba (t,Q =e -«/» g^^^2) . 
A direct computation gives the following Duhamel type formula 

B(t,x) = C 1 (t)B°(x) + C 2 (t)(B°/2 + B 1 )(x) + [ £ 2 (t - s)curl (G)(s, x)ds, 

Jo 

with B 1 = d t B(t = 0) = —curl (E(t = 0)) = —curl (E°). As this was observed in 

there exists < c < 1 such that we have the following bounds 
• For |£| > 2 



l*i(t,OI < e" 
\$2(t,0\ < 



ct 

e ' 



For 1/4 < |£| < 2 

\$l(t,0\ + \$2(t,0\<e~ Ct 
For 2 9 " 1 < |f | < 2 q+1 , with q < -3 



|$i(t,£)l < :=e-2cosh(tVl/4-2%- 1 )) 
|$ 2 (t,OI < * 2 M):=e 



3 _ i sinh(^l-2 2 (^ 1 )) 
6 Vl/4- 22(9-1) ' 



On the one hand, for g > — 1, one has 

(5.3) \\A q B(t)\\ L 2 < e-^IIA^II^+e-^dlA^II^ + HA^ 1 !!^) 



t 

e- c{t - s) \\A q G\\ L 2ds. 







Multiplying both sides of (15. 3p by 2^2 applying Young's inequality (in time) 
and summing in q yields 



(5.4) \\(I-S )B\\ ,_ x < ||(I -S )5°||.$ _, + !!(/ -50)^11 ,_ a 

_L H H H 

+ \\(I-S )G\\ 2 a_ x 
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On the other hand, for q < 0, one has 

\\A q B(t)\\ L 2 < ^(OIIA^II^ + ^^^HA^II^ + IIA^ !!^) 

+ 2* [ $ 2 q (t-8)\\A q G(8)\\ L2 d8. 

Jo 

Taking the L 2 ^ norm in time and applying Young's inequality we get 

\\\B\\ L 2 L 2 < ||<I>^|| L 2( R +)|[A 9 S || L 2 

+ H^II^K+^HAgB 1 !!^ + \\A q B°\\ L 2^ 

+ 29 |l*glli 1 (R+)ll^qG ! |li|Z, 2 - 

But since for every q < and r € [1, +oo], & l q satisfies 

W g \\ L r(R + ) < 2~' q , i = 1,2. 
Multiplying both sides by 2 q 2q a / 2 and taking the I 2 norm gives 

(5.5) \\S B\\ A < \\S B°\\,a 1 + ||5 S 1 ||. # _ a + \\S Q G\\ 

Putting together (|5.4p and (|5.5j) gives (|3.3p as desired. □ 

Proof of Proposition 3.5 The proof is based on the paraproduct decompo- 
sition. We choose to prove in details only estimates f|3.6[) and (|3.7p . The other 
estimates are easier or classical and left to the reader. 

Proof of (|3.6p We decompose EB into 

EB = T E B + T B E + S 2 R(E, B) + (I - S 2 )R(E, B), 
and will show the following estimates: 

(5.6) \\T E B + T B E\\ i2T g-i {R2) < \\E\\ L 2 tL 2 {r2) \\B\\ L¥L2{r2) 

( 5 - 7 ) \\ S 2R(E,B)\\ l 1 tL 2 < \\E\\ L 2_ L2(R2) 11-811^2 ^n, 0(R 2) 

(5.8) - S 2 )R(E,B)\\ L 2 T ,g-i {R 2 ) < \\E\\Liq om {^?)\\ B \\L^q o& 

First, we prove (|5.6p . Since the term T B E can be treated in a very similar way, we 
focus on TeB. First, 

A q (T E B)= £ A q (A q BS q E). 

|5-«|<1 

Since A q is uniformly bounded on L 2 , we have 

^2-«||A,(r £ B)|| i4ia <X;2- a £ IIA^^II^. 

9 9 |g-g|<l 

We are going to deal with the term q = q only (the two other terms q = q db 1 can 
be estimated similarly). Applying successively Holder's inequality (in the variables 
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t and x), Bernstein's lemma, Young's inequality (in the variable q), and Holder's 
inequality (in the variable q) gives 

Y,^ q \\^ q BS q E\\ LlL , < ^2-«||A 9 S||^ i2 ||S 3 £||^ 

^ E 2 ~ 9 E 11 A ^n^ 2 ii ajew^ 

i j<q 



<EE 2J " 9 ii A ^ii^ 2 ii A ^n^ 2 




Next we prove (|5.7p . Applying Bernstein's Lemma [37X1 and Cauchy-Schwarz (in j) 
gives 

||S 2 5(5,5)|| £l ^ 2 < ^11^(5,^)11^ 

<?<o 

< £>||A g 5(5,5)||^ Ll 

<J<0 

< j> E iia^ii^iia^ii^ 

^ E E ^iia^ii^iia^ii^^ 

j g<inf(0,j+2) 

< ^2^|A i 5|| L|i2 ||A i 5|| i|i2 +^||A i 5|| i|L2 ||A i 5|| I , |i2 
i<o i>o 

< \\e\\l* t l* II^II^lo- 

To estimate (|5.8j) . Holder's inequality (in t, x) and Cauchy-Schwarz (in j) gives 

||(J-5 2 )5(5,5)|| £| g_i < E E IIA^II^HA^IUcc^ 

g>0j>g-2 

^ E E \\*jE\\i%i*\\AjB\\l?& 

j>~20<q<j+2 

< Yl maxO^^IIA^II^IIA^Hioo^ 
i>-2 

< \\E\\L*q \\B\\ L ™q ■ 

T log T log 

Proof of (|3.7p As for the proof of (|3.6p , we split uB following the paraproduct de- 
composition: 

uB = T B u + T U B + R(u, B). 
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We shall only estimate here Tbu, the estimate of R(u, B) being similar, and that of 
T U B easier. By Holder's inequality, 



\ T Bu\\h T Ll g =X)max(l,g) \\S q BA q uf L 2 L 2 
q 

< ^max(l,g) \\A q u\\ 2 L 2 L 2 \\S q B\ 



Now observe that Bernstein's lemma and Cauchy-Schwarz' inequality (in j) give 

\ 1/2 



maxfl, q) 



< 2Q 



v /max(l,g) 

Coming back to the bound for ||Teii|| L 2 L 2 , this gives 

T log 

ii T ^iii^ og <E 229 ii A ^ii^ii^iiW 



References 



[i 
[2 

[3 
[J 

[5] 
[6] 

[7 

[8 

[9 

[10 

[11 



[12 
[13 



H. Bahouri, J.-Y. Chemin, and R. Danchin. Fourier Analysis and Nonlinear Partial Differential 
Equations, Grundlehren Math. Wiss. 343, Springer- Verlag, Heidelberg, 2011. 
D. Biskamp. Nonlinear magnetohydrodynamics, volume 1 of Cambridge Monographs on Plasma 
Physics. Cambridge University Press, Cambridge, 1993. 

J.-M. Bony. Calcul symbolique et propagation des singularites pour les equations aux derivees 

partielles non lineaires. Ann. Sci. Ecole Norm. Sup. (4), 14(2):209-246, 1981. 

M. Cannone, Q. Chen, and C. Miao. A losing estimate for the ideal MHD equations with 

application to blow-up criterion. SIAM J. Math. Anal, 38(6):1847-1859 (electronic), 2007. 

J. -Y. Chemin, Perfect incompressible Fluids, Clarendon press, Oxford, 1998. 

J. -Y. Chemin and I. Gallagher. On the global wellposedness of the 3-D Navier-Stokes equations 

with large initial data. Ann. Set. cole Norm. Sup. 39(4), no. 4, 679698, 2006. 

J.-Y. Chemin and N. Lerner. Flot de champs de vecteurs non lipschitziens et equations de 

Navier-Stokes. J. Differential Equations, 121(2):314-328, 1995. 

P. A. Davidson. An introduction to magnetohydrodynamics. Cambridge Texts in Applied Math- 
ematics. Cambridge University Press, Cambridge, 2001. 

H. Fujita and T. Kato. On the Navier-Stokes initial value problem I, Archiv for Rationnal 
Mechanic Analysis, 16, (1964), pages 269315. 

S. Ibrahim and S. Keraani. Global small solutions of the Navier-Stokes-Maxwell equations, To 
appear in SIAM jour. Math. Anal.. 

S. Ibrahim and T. Yoneda. Local solvability and loss of smoothness of the Navier-Stokes- 
Maxwell equations with large initial data , http://arxiv.org/pdf/1109.6089.pdf. To appear in 
JMAA. 

P. G. Lemarie-Rieusset. Recent developments in the Navier-stokes problem. CRC Press, (2002). 
J. Leray. Essai sur le mouvement dun liquide visqueux emplissant lespace, Acta Mathematica, 
63, (1933), pages 193248. 



16 



PIERRE GERMAIN, SLIM IBRAHIM, AND NADER MASMOUDI 



[14] N. Masmoudi. Global well posedness for the Maxwell-Navier-Stokes system in 2D. Journal de 
Math. Pures et Appl. 93 (2010) 559571. 

The Courant Institute,, New- York University, Mercer Street, USA 
E-mail address: pgermain@cims.nyu.edu 

Department of Mathematics and Statistics,, University of Victoria, PO Box 3060 
STN CSC, Victoria, BC, V8P 5C3, Canada 
E-mail address: ibrahim@math.uvic . ca 
URL: http : //www. math. uvic . ca/~ibrahim/ 

The Courant Institute,, New- York University, Mercer Street, USA 
E-mail address: masmoudi@cims.nyu.edu 



